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UVOT Measurements of Dust and Star Formation in the SMC and M33 


1. Introduction 

Finding the rate at which galaxies form stars is fundamental to understanding the galaxies’ 
formation and evolution. One can measure the star formation rate (SFR) at many wavelengths, 
using both broadband and line-based methods (see, e.g., Kennicutt & Evans 2012). Each of these, 
of course, have their strengths and weaknesses. 

Ultraviolet (UV) light traces the photospheres of O and B stars, and measures star formation on 
timescales of ^100 Myr. However, it is more heavily extinguished by dust than longer wavelength 
indicators, so it is necessary to do an accurate dust correction. Unfortunately, the shape of the 
dust extinction curve is uncertain in UV. Figure 1 shows four commonly used extinction curves, 
which are nearly identical at optical wavelengths, but vary widely in the UV. There are two main 
differences in the extinction curves: the slope (typically parametrized as Ry) and the strength of 
the bump at 2175 A. 



Figure 1: Several common dust extinction curves: Milky Way (Cardelli et al. 1989), starburst (Calzetti et al. 
2000), Large Magellanic Cloud (LMC; Misselt et al. 1999), and Small Magellanic Cloud (SMC; Gordon 
et al. 2003). Note the extinction curves’ differing slopes and the presence (or lack) of the 2175 A dust 
bump. Overlaid are the filter transmission curves for the uvw2, uvm2, and uvwl filters, showing how their 
placement is useful for constraining the properties of the dust extinction curve. 

Previous measurements of the extinction curve have relied on the time-consuming procedure 
of observing the spectrum of one extinguished star and comparing it the spectrum of an unextin¬ 
guished star of identical spectral type. It is unsurprising, therefore, that very few extinction curve 
measurements have been done in this way. The standard curve associated with the Small Magel¬ 
lanic Cloud (SMC), for example, is based on five lines of sight (Gordon et al. 2003), and the other 
curves rely on similarly few measurements. 

Our goal is to model the dust extinction curves in the SMC and M33 and determine if the 
curves vary across the galaxies. The UV/optical telescope (UVOT; Roming et al. 2005) on Swift 
(Gehrels et al. 2004) has three near-UV filters, which are included in Figure 1. The uvm2 filter 
overlaps the 2175 A dust bump, and the uvwl and uvw2 filters are to either side, which means 
we can constrain the dust curve. This is demonstrated further in Figure 2: comparing the unex- 


2 














UVOT Measurements of Dust and Star Formation in the SMC and M33 



1500 2000 2500 3000 3500 4000 4500 
Wavelength 


Figure 2: Model spectra representing the ability of UVOT to measure dust extinction properties. The top 
spectrum has no dust attenuation (Ay = 0). The bottom spectrum has a small attenuation of Ay = 0.15 with 
a dust extinction curve similar to that of the Milky Way. The UVOT fluxes for each spectrum are marked 
with red diamonds. Compared to the top spectrum, the uvm2 flux is suppressed relative to uvwl and uvw2 
(which traces the 2175 A bump), and the uvwl and uvw2 fluxes are more similar (which traces Ry). 

tinguished spectrum to the extinguished one, the uvm2 magnitude is suppressed relative to uvwl 
and uvw2, and a line connecting uvwl and uvw2 becomes shallower. The degree to which these 
happen corresponds to the strength of the 2175 A bump and the value of Ry, respectively. 

To accomplish this goal, we use a mosaic of UVOT and ground-based optical observations of 
the SMC and M33, the details of which are described in Section 2. We model the spectral energy 
distributions of star-forming regions in the galaxies in Section 3, and describe our early results in 
Section 4. 

2. Data 

We observed the SMC with UVOT as part of the SUMaC (Swift Ultraviolet survey of the 
Magellanic Clouds) program. The SMC observations consist of 50 individual tiles, each 17 x 17 
arcminutes, and a total exposure time of 1.8 days. The M33 observing campaign consisted of 11 
hours of exposure divided between 13 tiles. False-color mosaics of the SMC and M33 UV imaging 
are in Figure 3. The UVOT data reduction followed the procedure discussed in Siegel et al. (2014). 

We used archival ground-based optical imaging of both the SMC and M33. The SMC optical 
data were from Massey (2002). The data were taken at the Curtis Schmidt telescope at CTIO in 
the Harris UBVR filters (Massey et al. 2000). Due to difficulties calibrating the [/-band data, we 
discarded it from analysis. The M33 data were from Massey et al. (2006), taken at the Mayall 
telescope at KPNO in the Harris UBVRI filters. 

All of the UVOT and optical images were aligned and rebinned to 10" (2.9 pc) pixels for the 
SMC, and to 2.45" (10 pc) for M33, using SWarp (version 2.19.1; Berlin et al. 2002). To remove 
the diffuse background in each image, we used a circular median filtering technique following 
Hoversten et al. (2011). Pleuss et al. (2000) used HST imaging of HII regions in MlOl to show 
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Figure 3: False color UVOT images of the SMC (left) and M33 (right). The image of the SMC (M33) is 
about 2.3° (1.2°) across, corresponding to 2.4 kpc (17.6 kpc) at the assumed distance of 60 kpc (840 kpc). 

that they range between 20 pc and 220 pc in size; we therefore adopt this range of values for the 
circles’ diameters. 

Finally, we used Source Extractor (SE; version 2.5.0; Bertin & Arnouts 1996) to identify star¬ 
forming regions in the UVOT uvw2 image. We chose this filter because it is the bluest available, 
thus tracing the most massive young stars. (The uvw2 filter has a red leak, but the transmission 

o 

isn’t significant beyond ^3000 A.) SE has a known problem in which it can identify pixels in 
non-contiguous regions as belonging to the same region; to correct for this, we follow the method 
described in Appendix A of Hoversten et al. (2011). After this correction, we find 396 star-forming 
regions in the SMC and 656 in M33. We then use the SE-defined regions to extract photometry in 
the other UV and optical filters. 

3. Modeling 

We model the spectral energy distributions of each star-forming region by comparing our data 
to a grid of models. We create these models using the PEGASE.2 spectral synthesis code (Eioc 
& Rocca-Volmerange 1997). The grid includes spectra for ages of 1 Myr to 13 Gyr, total dust 
attenuation (Ay) from 0 to 7 magnitudes, dust extinction curve slopes (Ry) from 1.5 to 5.0, and 
2175 A bump strength from 0 to 2 (where 1 is the strength found in the Milky Way). We assume 
that since we have isolated individual star-forming regions, a single starburst is a reasonable star 
formation history. 

At each point in the grid, we compared the colors of each star-forming region to those of the 
model and calculated a chi-squared. By comparing colors, and not magnitudes, we could account 
for the mass normalization at a later stage. We repeated this 250 times while moving the magnitudes 
within their uncertainties. We then created a distribution of the best fits from each of the 250 runs, 
from which we extracted an overall best fit and uncertainties. The mass normalization was derived 
based on this best fit. 

In upcoming papers, this chi-squared procedure will be replaced with a more statistically ro¬ 
bust Markov Chain Monte Carlo fitting technique. This will enable better measurements of degen¬ 
eracies and more reliable uncertainties. We will also fit the magnitudes, rather than the colors, and 
include the mass as a free parameter. 
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Figure 4: Top row: Maps of dust properties for star-forming regions in the SMC. We assume that we don’t 
constrain Ry or the bump strength in regions where Ay is less than 0.1. Bottom row: Map of the ages of 
star-forming regions, and the corresponding star formation history. 

4. Early Results 

The best-fit properties of the star-forming regions in the SMC and M33 are mapped in Fig¬ 
ures 4 and 5. Also in these figures are the star formation histories derived from the regions’ masses 
(not shown) and ages. 

The most important result is that the dust extinction curve varies significantly across both the 
SMC and M33. For studies of objects within these galaxies, the best-fit extinction curve can be 
used in place of others. In particular, the SMC curve (plotted in Figure 1) does not appear to 
apply to young star-forming regions in the SMC. For most other dust corrections in the literature, 
by necessity, authors choose their favorite extinction curve with little ado. This work shows how 
variable the extinction curve is, and therefore can be used to quantify the systematic error when 
correcting for dust extinction. 

The second important result is the detailed star formation history of the star-forming regions 
over the past few hundred million years. When using optical data, it is difficult to distinguish 
between the most massive stars, because it’s already in the Rayleigh-Jeans tail. As such, measure¬ 
ments of the youngest ages are difficult. By utilizing ultraviolet light, we can better resolve the 
ages of the youngest populations, and better understand the most recent star formation. 

In upcoming work, we will also be mapping the dust properties and ages across the entirety of 
both galaxies. Although the resolution will be coarser (due to computational constraints), this will 
yield useful information for more than just the brightest regions of star formation. 
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Figure 5: Same as Figure 4, but for M33. 
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